Vertex functions for d-wave mesons in the light-front approach 
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Abstract 

While the light-front quark model (LFQM) is employed to calculate hadronic transition matrix 
elements, the vertex functions must be pre-determined. In this work we derive the vertex functions 
for all d-wave states in this model. Especially, since both of 3 Z?i and 3 «Si are 1 mesons, the 
Lorentz structures of their vertex functions are the same. Thus when one needs to study the 
processes where 3 D\ is involved, all the corresponding formulas for 3 S*i states can be directly 
applied, only the coefficient of the vertex function should be replaced by that for 3 D\. The results 
would be useful for studying the newly observed resonances which are supposed to be d-wave mesons 
and furthermore the possible 2S-1D mixing in ip' with the LFQM. 
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I. INTRODUCTION 



In the hadronic physics the most tough job is to calculate the hadronic transition ma- 
trix elements which are fully governed by non-perturbative QCD. Because of lack of solid 
knowledge on non-perturbative QCD so far, one needs to invoke phenomenological models. 
Applications of such models to various processes have achieved relative successes so far. 
Among these models the light-front quark model(LFQM) is a relativistic model and has ob- 
vious advantages for dealing with the hadronic transitions where light hadrons are involved 



H, |2J. The light-front wave function is manifestly Lorentz invariant and expressed in terms of 
fractions of internal momenta of the constituents which are independent of the total hadron 
momentum. This approach has been applied to many processes and thoroughly discussed in 
literatures j3-13|. Generally the results qualitatively coincide with experimental observation 
and while taking the error ranges into account (both experimental and theoretical), they can 
be considered to quantitatively agree with data. 

However earlier researches with the LFQM only concern the s-wave and p-wave mesons 
whereas the higher orbital excited states have not been discussed yet. With the improve- 
ments of experimental facilities and rapid growth of database many new resonances have 
been observed and some of them are regarded as higher orbital excited states, for exam- 
ple ^(3770) and ^(4153) are suggested to be d-wave charmonia whose principal quantum 



numbers are respectively n — 1 and n = 2[l4f . New analysis on X(3872) indicates that 
it has a possible charmonium assignment 1 D 2 14J. Moreover, a d-wave state T(1 3 D 2 ) was 
observed by CLEO 15]. The situation persuades us to extend our scope to involve d-wave. 
It is generally believed that only the lattice theory indeed deals with the non-perturbative 
QCD effects from the first principle. So far, the lattice study is constrained by not only the 
computing abilities, but also the theory itself. Even so, remarkable progresses have been 
made on the two aspects. It is hoped that the lattice calculation will eventually solve all the 
problems on hadrons, such as the hadron spectra, wavefunctions, even the hadronic tran- 
sition matrix elements. At present, the lattice calculation indeed shed some light on the 
wavefunctions 
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For examples, the authors of Ref. [16[ suggested a method to compute 
the spectra and wavefunctions of hadron excited states and applied their method to the 
[/(l)2+i lattice gauge theory. Abada et. a/.[17| study the pion light-cone wave function on 
the lattice by considering the three-point Green functions. On other aspects, for example, 
the readers who are interested in the unquenched lattice calculation on the charmonia whose 
total quantum numbers J PC are determined simultaneously, are suggested to refer to those 
enlightening papers jl9| . Some works about the radiative decays of charmonia are included in 
concerned references [2 Ol. l2l| . More similar works can be found in most recent literatures 
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24J. In fact, due to the rapid progress of lattice calculations, people are more tempted to 



trust those results, but it by no means implies that we should abandon phenomenological 
models because those models are directly invented to manifest the physical mechanisms and 
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moreover, they are simpler and applicable in practice. 

To evaluate the transition rate in the LFQM one needs to know the wave functions of 
parent and daughter hadrons. For any 2S+1 Lj state, its wavefunction is constructed as the 
corresponding spinors multiplying the so-called vertex function which should be theoretically 
derived. It is also noted the wavefunctions for s-wave and p-wave have been derived and their 
explicit forms are given in Ref.[7|. Following the same strategy we obtain the wavefunction 
for all d-wave states. 

The traditional LFQM was employed to study the decay constants and form factors of 
weak decays but to maintain the Lorentz covariance other contributions such as Z- 

28| contributions must be included. 

6[ has been suggested which systematically includes the zero- 



25- 



diagram|5[ or zero-mode 

Thus a covariant LFQM 
mode contributions. In the traditional LFQM the constituent quarks in the bound state 
are required to be on their mass shells, nevertheless in the covariant LFQM approach the 
constituent quarks in the meson are off-shell while only the meson is on its mass shell. In this 
approach, one writes down the transition amplitudes where all quantities and the integrations 
maintain their four- dimensional forms, then integrates out the light-front momentum p~ in 
a proper way. While carrying out the contour integration the antiquark is enforced to be 
on its mass shell. The integrand in the remaining three-dimension integration reduces into 
a form where all quantities can be expressed in terms of the conventional wavefunctions. 
During this procedure, some extra contributions emerge comparing with the original scheme 
(see the text for details). 

In this work after this introduction we derive the phenomenological vertex functions for 
d-wave in the conventional light-front approach in section II. Then in section III we present 
their forms in the covariant light-front approach. In section IV we discuss some formula for 
3 Di states and the section V is devoted to a brief summary. 



II. VERTEX FUNCTIONS IN THE CONVENTIONAL LIGHT-FRONT AP- 
PROACH 

Let us first derive the vertex functions in the conventional light-front approach. 



In the conventiona 



J can be written as 
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light-front approach a meson with the total momentum P and spin 



\M(P, 2S+1 Lj, J z )) = J {d 3 p 1 }{d 3 p 2 }2(2TT) 3 5 3 (P - pi - p~ 2 ) 

E (Pi ) #2,Ai,A 2 )| gi (p 1 ,A 1 )g 2 (p 2 ,A 2 )), (1) 
A1A2 

where the flavor and color indices are omitted; q\ and g 2 correspond to the quark and 
antiquark in the meson and pi and p 2 are the on-shell light-front momenta of quark and 
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antiquark, p is the three-momentum (pi — p2)/2 and we define 

i 2-, - _m 2 + p 2 ± 3 _ _ dp+dV 



Pi = (pf,Pi±), Pix = (Pi,Pi), Pi = — -^,{ d Pi}- 2 (2tt) 3 ' 
l?i(pi, Ai)g 2 (P2, A 2 )) = &I>i)4 2 (p 2 )|0>. 
The light-front momenta px an d £>2 are expressed via the variables p and a?j (i = 1, 2) as 

p+ = XiP + , P2 = X 2 P + , Xi + X 2 = 1, 

Pi± = x 1 P ± + p ± , p 2 ± = x 2 P± - p±. (2) 




FIG. 1: The vertex for meson-quark- antiquark. 



In the momentum representation, the wavefunction ^ J L J g for the state 2S+1 Lj can be 
decomposed into the form 

^ J L f{p~ 1 ,p~ 2 ,X 1 ,X 2 ) = -^=(LS;L z S z \LS;JJ z )R s x ^(x,p ± )^ LLz (x,P±), (3) 



where <Pll z ( x iP±) describes the relative momentum distribution of the quark (antiquark) 
in the meson and L is the orbital angular momentum between the constituents. R\^\ 2 
transforms a light-front helicity (Ai,A2) eigenstate to a state with definite spin (S,S Z ) and 
it is expressed as 

< S a 2 (^) = T^TJV} \ , MpiAiW + M )T s v(p 2 ,X 2 ), 

V2M (M + mi + m 2 ) 

T = 75 (for£ = 0), 

ri =-${S z ) (forS = l), (4) 

with P = pi + p 2 , Mq = m ' 1 ^ ± + m ' 2 ^ ± and M = yj Mq — (mi — m 2 ) 2 . For more details, 
readers are suggested to refer Ref.P 
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In the LFQM, the harmonic oscillator wavefunctions are employed to describe the relative 
3-momentum distribution of quark and antiquark in a meson. For d-wave the harmonic 
oscillator wavef unction is 29] 



(5) 



where K = (p 2 —pi)/2 and ip is the harmonic oscillator wavefunction for s wave and its 
explicit expression is with 



<P = 4 (^) S/4 1 



dPz 
dx 2 

x 2 M ml + p\ 



Vz 



p 2 z + pI - 

2/3 2 ' 



2 2x 2 M 

Substituting Eq.(j3]) into Eq.flS]), we deduce the expression 

1 



(6) 



*2S (Pl,#2,Ai,A 2 ) 



(2S;L z S z \2S;JJ z )R s x X(x, PA _)cp 2Lz ( x , PA _) 



r N~ c V2M (M + m 1 +m 2 



-u(p u Ai)(f + M )T { 2s +1Dj) v( P2 , A 2 ),(7) 



with 



r 



(3^) 



CD2) 



( 3 D 2 ) 



( 3 D 3 ) 



5(3 4 



Ml 



Mq 



3/3 



A (K ■ P) 2 — M 2 K 2 



9/3' 



3M 

-i^J^K^g^ + KpK a g vP + K a K u g^) 



where relations (2 0; L z 0|2 0; 2 J z ) 



£*^(L,)^(J Z ), {21;L Z S Z \21 ] U Z ) 



*-e*^{L z )E^{S z )e v {J z ) , (2 1; L z S z \2 1; 2 J 2 ) 



and (21;L,£ 2 |21;3J Z ) = §£^V*)[£V(^K*<*(^) + eVW^^) + 
are used. 

One can further simplify these wavefunctions in terms of the Dirac equation fiu{jp\) = 
miu(pi) and pV^ (pi) = —m 2 v(pi), so that all scalar products of vectors are replaced by only 
Mo, mi and m 2 via a simple algebra, thus the wave function is 



^2S Z (P^P2, Ai, A 2 ) = u(pi, \ 1 )h' { 2s+i Dj) T' { 2s+i Dj) v(p 2 , X 2 ) 



(9) 
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where 





(1Z)2) \j N C M ^ 

= V^^' (10) 



and 



with 



W(3 Dl ) 

r' (lDa) = ^k^kj^, 

W ( a D 2 ) W eD 2 ) W { Z D 2 ) W ( 3 D 2 ) 

2K 2K 2K 

r' (3D3) = [# M jr„( 7a + — 2L) + i^ a (7, + — -) + KaKyfrp + — (n; 

1 3> W { 3 D:j) W { 3 D:j) W { 3 D:j) 



W(3 Dl) 



(mi + m 2 ) 2 - Ml 
2Mq + mi + m 2 



W (*D 2 ) 

w b {3D: 



-12M 2 



[M 2 - (mi + m 2 ) 2 ] [M 2 - (mi - m 2 ) 2 ] 

6M| 

(2M + mi + m 2 )[M 2 - (mi - m 2 ) 2 ] 



W(3 



6M 2 



W/3 



(3Da) (M - mi - m 2 ) [M 2 - (m x - m 2 ) 2 ] ' 
M 



(3D 2 ) ~~ ' 

m 2 — mi 



W(3 D 3) = M + mi + m 2 . (12) 



It is interesting to ask where the QCD which definitely governs the physical processes, gets 
involved or how can one implant the QCD information into our calculation in the LFQM. 
In Ref. 30J the authors derived an effective Hamiltonian for bound states in the light-front 
frame based on the standard Lagrangian of QCD. The vertex function is the effective coupling 
between the bound state and its constituent quarks, thus as the wavefunction of the bound 
state is obtained the effective vertex function is in hand. After a long discussion about the 
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-1>> 



FIG. 2: The Feynman diagram for a meson annihilation. 
Lorentz structure and the features of the dynamics of the vertex functions, they derive the 



vertex function which has exactly the form of Eq.Q in our work. In Ref . [30| . the radial wave 
function was obtained by solving the eigenvalue equation numerically. Obviously all QCD 
information (both short- distance and long-distance effects) is involved in the Hamiltonian 
and as well as in the solution. As argued in literature [5|, |29|, the solution can be well 
approximated by a Gaussian function with model parameters to be fixed by fitting data. 
Thus following Ref.0] we choose a Gaussian wave function where the QCD information is 
included in the model parameter (3. 

We can apply these wavefunctions to deal with concrete physical processes. For example, 
when we calculate the rate of 3 D± state annihilation through a vector current (in Figj2]), the 
transition amplitude is written as 



= N c I **f> "gL TYfo - ~ m,hM + m.)K- (13) 

III. VERTEX FUNCTIONS IN THE CO VARIANT LIGHT-FRONT APPROACH 

Comparing with the conventional LFQM where both pi and P2 are on their mass shells, 
in the covariant light-front approach the quark and antiquark are off-shell, but the total 
momentum P = p\ + p 2 is the on-shell momentum of the meson, i.e. P 2 = M 2 where M is 
the mass of the meson. Obviously, the covariant LFQM is closer to the physical reality. 

If one tries to obtain the covariant vertex functions based on an underlying principle, i.e. 
QCD, he should invoke a reasonable theoretical framework. To directly obtain the covariant 
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rors of Ref . |6| suggested to 
, [32)]. The kernel of the 



vertex functions in the 4-dimensional momentum space, the aut 
solve the Bethe-Salpeter (B-S) equations for the bound states 
B-S equation includes the Coulomb piece which is induced by the one-gluon exchange as 
well as its higher-order corrections, and the confinement piece which incorporates the non- 
perturbative QCD but is not derivable so far. Genrally for solving the B-S equation, the 
instantaneous approximation is usually taken. 

In the concrete calculations of the observable physical quantity in terms of the LFQM the 
final result is eventually reduced into an integration over the four-momentum. Fortunately, 
by doing so, we may not really need the explicit covariant wavefunctions defined in the 
four-momentum space. Namely, we try to reduce the integration into a simple form where 
only three-momentum wavefunctions remain by a mathematical manipulation, then we are 
able to relate the corresponding integrand to the conventional vertex function which is well 
defined in the three-momentum space. 

Since the Lorentz structures of the covariant vertices are the same as that of the conven- 
tional vertex functions we rewrite these covariant vertex functions in Eq.(j9]) as 

iH { s Dl) [^ - — ^ — (pi - j^e", 
iH { , D2)1 ^K v e^, 

iRpihM^—'f'f + 77^— 1^ + 777^7^ + 777^—^1^, 

VV ( 3 D 2 ) VV ( 3 D 2 ) VV ( 3 D 2 ) VV ( 3 D 2 ) 

IK IK IK 

iH ( s Ds) [K,K„( la + —2-) + K,K a { lv + — + K a K v {l, + 777^)]^ (14) 

W(3£> 3 ) W(3D 3 ) W( 3 As) 

where H(2s+i Dj \ and Wns+i Dj -\ are functions in the 4-dimensional space. Practically, the 
vertex function(s) is(are) included in a transition matrix element, for example, the amplitude 
of 3 Di state annihilation via a vector current is written as 

Acov -2 N c f ,4 -^( 3 Oi)rp rr (Pl-P2)u u / . \ / / . u v 



16tt 4 J ^ N±N 2 

where s = Tr{[y v — ^ r ^](-/ 2 + wi2)7/i(/i + m i)}, Ni = p\ — m\ + ie and N 2 = p^—ml + ie. 



One first needs to integrate over p 1 as discussed in Ref . |6|, |7J . Integrating over p 1 is completed 
by a contour integration where the antiquark is set on shell. Then the integration turns into 

N r r . .0 hi 



h / w57 ^, (16) 



16ir 3 J 1 r - L x 2 x 1 (M 2 - M ( 



where W(3 Dl ^ and i v replace W^ Dl ) and e v in Eq.( |T5i) respectively. 
Following Ref. [3] we have the relation 



h { 3 Dl) = (M 2 - Ml)^T 2 h' (ADi) . (17) 
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An additional factor (M 2 — Mq )sJx\X2 was introduced when comparing the decay constant 
fp of pseudoscalar meson obtained in the two approaches as depicted in the appendix A 
of Ref.0]. The legitimacy is guaranteed because the decay constant is free of zero mode 
contribution. Then the authors have applied the relation into the vertex functions for S and 
P waves. To show the reasonability of such replacement, we substitute Eq. (fTT|) into Eq. (fT6l) 
to obtain a new expression whose form is similar to the right side of Eq. ffT?]) . However, 
the trace in Eq. (TIB]) involves the zero mode contribution which makes its form different 
from that in Eq. lfiUI) . Generally, after the contour integration over p{ , h(2s+i Dj \, w^s+i Dj ^ 
and e replace H^s+i Dj ), W^s+i Dj ^ and e respectively with the following relation to the 
corresponding quantities of the conventional LFQM 

h { 2 S+ i Dj) = (M 2 - Ml)^x~ 2 h' ( 2s + i Dj) . (18) 

Here we only concern the form of the covariant vertex function for the D-wave, includ- 
ing its Lorentz structure and coefficient, as well as its relations to the conventional vertex 
function. The details about the S- and P wave vertex functions were discussed in earlier 
literature 0]. When one needs to calculate a transition matrix in the covariant light-front 
quark model, he must know those vertex functions. Jaus has analyzed the case of the covari- 
ance of the transition matrix jsj, and in his work, a general form of vertex function is used 
and the three- momentum conservation is automatically guaranteed. He j6| indicates that 
the general form of the vertex function h must be functions of = Xi(M 2 — Mq) (i = 1,2). 
Obviously the function h adopted in this work coincides with this requirement. 



IV. THE FORMULA FOR 3 D 1 STATE 



For a J PC = 1 state, the the orbital momentum between the two constituents may be 
L = (s-wave) or L = 2 (d-wave) and their total spin is 1 (S — 1). In Ref. 0] the authors 
gave the meson-quark-antiquark vertex for 3 Si state as 

iHvbt-^-fa-pzU (19) 
Carrying out the contour integration over p{, H v and Wy turn into hy and uy 

Wy = M + mi + "^2, 

where the subscript V only refers to 3 S\ state. 

The Lorentz structure of the vertex functions for 3 Di and 3 S\ states are the same because 
they have the same quantum number J PC . The difference between the s-wave and d-wave 
is included in the coefficient functions Km and %. 
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The decay amplitude of an Si state via a vector current is proportional to 

A, = -i 2 ^ J d'pi^-Tr^ + mi )[ lv - ij ^ k ](-fc + m*)}?, (20) 

which is the same as Eg. ( TT5|) . except Hs Dl and W3 Dl are replaced by Hy and Wy respectively 
Integrating over p\ Hz Dl , Wa Dl , H v and Wy reduce into h3 Dl , ws Dl , h v and wy. The decay 
constant for 3 Si state reads 

Iv = ^mS dX ^ XlX2 {M* ~ Mg) [XlM ° 2 - mi(mi - m2) - Pl + "^W^^ (21) 

so that one would obtain the decay constant of the 3 Di state by replacing hy and wy by 
h-i Dl and ws Dl , thus it is 

f N c f J J2 ha D! r A/ r2 ( \ 2 . ™l+™2 2 1 /oo\ 

/3Bl = 4^M 7 ^(MS-Mg) ^ - mi(mi " m2) + ^T Pj ' (22) 

In fact since the Lorentz structure of the vertex functions for 3 Di and 3 S± are the same 
all the formula for 3 D\ can be deduced from those for 3 Si. For example, the form factors 
f,g,a+ and a_ of P — > 3 Di(V) decay can be obtained by simply replacing hy and wy of 
f,g,a + and a_ given in Ref. a, |7|. 

With these formula we will be able to explore some new resonances of angular excited 
states, or furthermore to study the mixing of 25" — ID which was proposed to explain the 
famous p — Ti puzzle for ?//|33l-l35l] in this model. 



V. A BRIEF SUMMARY 

In this paper we deduce the vertex functions (or wave functions) for the d-wave in the 
conventional and covariant light-front quark model. 

For the 3 D\ state the J PC is 1 and the Lorentz structure of its wave function is the 
same as that for the 3 S'i state so we obtain some useful formula for 3 D\ from the formula 
for 3 S'i given in Ref. 0,0]. It is noted we just discuss the vertex functions (iT m) for the 
incoming meson whereas for the outgoing meson the corresponding vertex functions should 
be i(7oIm7o)0- 

Since we adopt the Gaussian-type function for the radial part of the whole wavefunction 
instead of a solution obtained by solving the Schrodinger equation or the B-S equation, the 
simplification definitely brings up certain theoretical uncertainties, but as more data will be 
collected in the future, the more precise model parameter(s) will be determined and even 
the form of the wavefunction can be improved, thus we may do a better job along the line. 

These vertex functions can be employed when one calculates the transition rates in this 
model. In the future we will study some concrete physical transitions where d-wave mesons 
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are involved in terms of these vertex functions. The results will be compared with data 
and the consistency would tell us the accuracy degree of the model and the derived vertex 
functions. Once the validity of the model is verified via some processes, we can further discuss 
some long-standing puzzles and help to identify new resonances which are continuously 
observed at BES and BELLE and elsewhere. 
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